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Abstract 

We consider many-body quantum systems that exhibit quantum chaos, 
in the sense that the observables of interest act on energy eigenstates 
hke banded random matrices. We study the time-dependent expectation 
values of these observables, assuming that the system is in a definite (but 
arbitrary) pure quantum state. We induce a probability distribution 
for the expectation values by treating the zero of time as a uniformly 
distributed random variable. We show explicitly that if an observable 
has a nonequilibrium expectation value at some particular moment, then 
it is overwhelmingly likely to move towards equilibrium, both forwards 
and backwards in time. For deviations from equilibrium that are not 
much larger than a typical quantum or thermal fluctuation, we find 
that the time dependence of the move towards equilibrium is given by 
the Kubo correlation function, in agreement with Onsager's postulate. 
These results are independent of the details of the system's quantum 
state. 
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I. INTRODUCTION 



Many-body systems typically exhibit certain dynamical properties that are studied under 
the subject headings of thermodynamics and statistical mechanics. These properties include 
the following: 

1) Given an arbitrary initial state, the system almost always evolves towards an identifi- 
able condition known as thermal equilibrium, and then remains in this condition at almost 
all subsequent times. 

2) When the system is in thermal equilibrium, observables of interest take on values that 
depend only on the nature of the system and its total energy, but not on any other details 
of the specific state of the system. 

3) When the system is in thermal equilibrium, the measured value of an observable of 
interest at any particular time fiuctuates about its equilibrium value, with fiuctuations whose 
amplitude is suppressed by a factor of A^~^/^, where N is the number of degrees of freedom. 

4) During the approach to thermal equilibrium, the values of observables of interest 
are governed by equations that are not time reversal invariant. These equations typically 
depend on the values of other observables, possibly at different times. The information about 
the system that is included in these equations is not sufficient to reconstruct the complete 
physical state of the system. 

5) Often (but not always), these equations are markovian; that is, they depend only 
on the values of the observables in question, and their first time derivatives, at any given 
moment. 

There is a vast literature on the derivation of these properties from an underlying de- 
terministic, time reversal invariant dynamics, classical or quantum. In this paper (closely 
related earlier work includes [|l|-|T^), we explore to what extent these properties can be de- 
duced as consequences of quantum chaos. This means that we will assume that the energy 
eigenvalues and (more importantly) the matrix elements of typical observables have certain 
properties. These properties are believed (and, in some cases, rigorously proven) to hold for 
a canonically quantized system whose classical phase space is fully chaotic at the energies of 
interest, and they are likely to hold at least approximately for a broader array of systems. 

The outhne of this paper is as follows. In Section 2 we state our basic assumptions, and 
briefiy discuss their origins in quantum chaos theory. Some previous work that is directly 
relevant is summarized in Section 3. Sections 4 and 5 present new results concerning the 
approach to thermal equilibrium. Section 6 discusses the main conclusions. 



II. QUANTUM CHAOS 

We assume that the quantum system of interest is bounded and isolated, with N degrees 
of freedom, where 1. Since the system is bounded, the energy eigenvalues are discrete, 
and since it is isolated, its time evolution is governed by the Schrodinger equation. Let 
denote the energy eigenvalue corresponding to the energy eigenstate the state of the 
system at any time t is then given by 

|V't) = Ec«e-^^"*|a). (2.1) 
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We emphasize that |a) is an eigenstate of the full, many-body hamiltonian. The c^'s specify 
the state at any one time (say, t = 0), and we assume the usual normalization 

{A\A)=J2\^a\' = ^- (2-2) 
a 

Note that we have set ^ = 1 to simplify the notation. We will, however, point out how 
various quantities scale with h when this information is useful. 
We now make two key assumptions about the system. 

Our first assumption is not strictly necessary, but it will simplify some of the subsequent 
analysis. We assume that if any two sums of equal numbers of energy eigenvalues are equal, 

+ . . . + Ea„ = Ef3^ + . . . + Ep^ , (2.3) 

then {Pi, . . . , Pn} is a permutation of {ai, . . . that is, the corresponding eigenstates 

are the same for both sides. In particular, E^ = Ep implies that |a) = so that all the 
energy eigenvalues are nondegenerate. This assumption is expected to hold in general for 
a quantized chaotic system in which all unitary symmetries (such as invariance under any 
discrete or continuous reflection or rotation) have been removed by suitable changes (such 
as by putting the system in an irregularly shaped box). In this case, the energy eigenvalues 
have the same statistical properties as the eigenvalues of gaussian random matrices |]T^ . The 



system may or may not be invariant under the anti-unitary transformation of time reversal; 
we will assume that it is, because one of the most interesting aspects of thermodynamics is 
the appearance of an "arrow of time" even when the underlying dynamics is time reversal 
invariant. 

Our second assumption is the crucial one. Let A be a hermitian operator corresponding to 
an observable of interest that is a smooth, ^-independent function of the classical coordinates 
and momenta. We assume that the matrix elements of A in the energy eigenstate basis take 
the form ||18|,|6| 



A^p = A{E)5^p + e-^(^)/V(^, i^)Rap . (2.4) 

All the factors in this formula need explanation. 
First, for notational convenience we have defined 

E = \{E^ + Ep) and uj = E^-E^. (2.5) 

S{E) is the thermodynamic entropy at energy E, given by 

e^(^) = E5:5,(E-E„), (2.6) 

a 

where 5e{x) is a Dirac delta function that has been smeared just enough to render S{E) 
monotonic. A{E) and f{E,uo) are smooth functions of their arguments whose physical 
implications will be the main focus of this paper. Finally, Rap is a numerical factor that 
varies erratically with a and 13. It is helpful to think of the real and imaginary parts of Rai3 
as random variables, each with zero mean and unit variance. Without loss of generality, 
we can take f{E,uj) to be real, positive, and an even function of uj; then hermiticity of A 
implies R^a = R*afi- Also, in many cases of physical interest, R^js is purely real. 
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Eq. ( |2.4|) is semiclassical in nature; the factor of e'^^^^^"^ scales like h^^~^^^'^. Thus the 
validity of eq. ( ^.41 ) requires h to be "small," which in practice means that the energy E must 
be "large." The appropriate definitions of "small" and "large" are a key problem of quantum 



chaos theory. For systems with few degrees of freedom, it is now well established |^ that a 
necessary condition for quantum chaos is 5 <^ ^/Txh? where 6 ~ e~^E is the mean spacing 
between energy eigenvalues near E, and rxh is the Thouless time [^] (roughly speaking, 
the time scale by which all diffusive classical processes have saturated). For many-body 
systems in general, less is known. Studies of nonlinearly coupled oscillators and of 



fermions with pseudo-random single-particle energies and two-body matrix elements |TT]|T^ 
both suggest a threshold energy for quantum chaos that goes to zero in the thermodynamic 
limit like N~'^, with < < 1. We may therefore hope that eq. ( |2.4| ) will enjoy a wide 
range of validity. 

An important feature of eq. ( p. 41) is that the general structure that it describes is pre- 
served under multiplication 0. Thus, for example, the matrix elements of any power of A 
are given by 

{A^U = An{E)6^p + e-'^^^/'UE, c.)/?^ , (2.7) 

where An{E), fn{E, uj), and R^^^ can be expressed in terms of A{E), f{E, uj), and Rap- The 
precise relationship will not be needed, however; the key point is the generic character of 
eq. (0). 

Finally, the function A{E) can be related to a standard expression in statistical mechan- 
ics: the equilibrium value of A, as given by the canonical thermal average 

Wt^^^^. (2.8) 

Here T is the temperature, and we have set Boltzmann's constant to one. To see the relation 
between {A)t and A{E), we use eqs. ( p^) and ( p.6|) in eq. ( p.8|) to get 



TOO dE S{E)-E/T AfTp\ 

(^}T^ '"g%^,,,,_^r +0(e-^'''). (2.9) 

Jo E 

When N is large, the entropy is extensive: S{E,N) = Ns{E/N) + O(logA^). In this case 
the integrals in eq. (|2.9|) can be evaluated by steepest descent. We see that their ratio 
is A{E) + 0{N~^), where E is now fixed in terms of T by the steepest-descent condition 
dS/dE = 1/T; this also implies E = {H)t- Turning around eq. ( p.9|) then gives us the 
expression for A{E) that we want 

A{E) = {A)t + O(iV-i) + 0(e-^/2) _ (2.10) 

We will always assume that is large enough to make the indicated corrections negligible. 
Note also that eq. ( |2.1CI| ) is consistent with, but not identical to, Shnirelman's theorem 



21 1 . This theorem essentially states that A{E) is given by the classical, microcanonical 



average of A, up to corrections which are expected to be 0(/i^^^). Eq. ( |2.10D , on the other 



hand, already includes corrections up to 0{h^^ ^''^^)j but has in addition corrections of 
0{N-'). 
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III. THERMAL EQUILIBRIUM 



The expectation value of an observable A in the state specified by eq. (|2.1|) is given by 



A = {tPt\A\iJt) 



Cq,C^ e /la/3 

a/3 



(3.1) 



We will take At as the main object of study. It is not obvious that this is the right thing to do, 
since short-time measurements do not generally yield quantum expectation values. However, 
our main goal is to compare with the results of conventional nonequilibrium statistical 
mechanics, in which time-dependent expectation values are the basic ingredients (see, e.g., 
P^^4|, and Section V, below). A detailed discussion of the quantum measurement problem 
would be needed to address this issue properly, but this is beyond the scope of the present 
paper. 

If we now take the infinite time average of At, we find 



A = lim 



r JO 



dt At 



E 



r P4 
'^a I ^aa ) 



(3.2) 



where the last line requires nondegeneracy of the energy eigenvalues. It also requires that 
the averaging time r be much larger than the Heisenberg time th = 27rh/6 ~ e"^. This time 
scale is much too large to be physically relevant, and thus the infinite time average must 
be regarded as a purely theoretical device. Nevertheless, if the system comes to thermal 
equilibrium, then At should be near its equilibrium value (^4)^ most of the time, and thus 
we should have A = {A)t- 

To check whether or not this is the case, we first substitute eq. (|2.4|) into eq. (p.2D to get 



A = J2\ca\'AiE^) + 0ie 



(3.3) 



We now make a mild assumption about the state \ipt)- We assume that the expectation 
value of the total energy 



(3.4) 



has a quantum uncertainty 



A 



nl/2 



(3.5) 



that is small, in the sense that A'^\A"{E)/A{E)\ ^ 1. This is a natural assumption when 
is large, since states of physical interest typically have A ~ N~^^'^E. If we now expand 
A{Ea) in eq. ( ^.3] ) in powers of E^ — E, we get A = A{E) + O(A^); combining this with 
eq. (|2.10|) , we find 
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A = {A)t + 0(A2) + O(Ar-i) ^ o(e-s/2) _ ^3_g) 

Thus we have shown that the infinite time average of At is indeed equal to its equihbrium 
value at the appropriate temperature. Note that this key property follows entirely from 
the matrix element structure of eq. (|2.4| ), and does not depend on the details of the initial 
quantum state. 

We now turn to an examination of the fiuctuations of At about its equilibrium value A. 
The mean squared amplitude of the these fiuctuations is 



{At -Af=\imJ- dt {At -Af 



T JO 

2 \„ \2\ A |2 



= 0{e~') . (3.7) 

We see that the fiuctuations of At about A are very small. This tells us that, whatever the 
initial value Aq happens to be, At must eventually approach its equilibrium value, and then 
remain near it most of the time. 

On the other hand, eq. ( |3.7| ) is too small to represent the expected thermal fiuctuations 
of A. To find them, we must look at what are usually called quantum fiuctuations. The 
mean squared amplitude of the quantum fiuctuations is 



{i^t\{A-AY\^l,t) = {A^)t-A^ 

= {A')t - {Af^ 

+ 0{A^) + 0{N-') + 0{e-^/^) . (3.8) 

In the last line, we have used the fact that the matrix elements of A"^ have the same general 



structure as the matrix elements of A, as shown in eq. ( p.7|) , and that this structure implies 
that the infinite time average is the same as the thermal average, as shown in eq. ( ^.61 ). 
Eq. ( p.8|) tells us that the quantum fiuctuations in A have the right magnitude to be identified 
as thermal fiuctuations 0. 

Note, however, that (A^)t — {A)t is itself expected to be 0{N'^^) for typical observables 
of interest (see, e.g., I^S])' ^^'-^ ^^^^ correction terms on the right-hand side of 



eq. (|3.8| ) are not necessarily smaller than the leading term. This is not a point of concern, 
however; we used the canonical ensemble to define thermal averaging, and the result would 
in general change by a factor of order one if we were to use instead the grand canonical or 
microcanonical ensemble. Since the exact size of the thermal fiuctuations in any particular 
observable depends on the choice of ensemble, our claimed identification of quantum fiuc- 
tuations as thermal fiuctuations can be meaningful only up to a numerical factor. This is 
what is established in eq. (p.8|). 

Another point of interest is the nature of the classical limit. Recalling that e""^ ~ 



we see that eq. ( |3.7| ) predicts that {At — AY vanishes in the classical limit. This is in 
fact a reasonable result if the classical system is chaotic. To see why, first note that the 
classical limit of a generic quantum state (to the extent that such a thing exists at all) 
is a probability density on phase space. Then the time dependent expectation value At 
is given (in the classical limit) hy At = J d'^^X Utp{X)A{X), where X is a point in phase 
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space, p[X) is the probability density associated with the quantum state at t = 0, and 
Ut is the Frobenius-Perron evolution operator for phase-space densities. When suitably 
regulated and renormalized, this formally unitary operator acquires eigenvalues (known as 
Ruelle resonances) inside the unit circle that are associated with irreversible decay to the 
ergodic distribution [^. If both A{X) and p{X) are smooth functions, At approaches a 
fixed limit as t — oo, and {At — AY vanishes. The mean squared amplitude of the classical 
thermal fluctuations is then given by the infinite time average of / cP^X Utp{X)[A{X) — A]^, 
which is the classical limit of the quantum expression {il)t\{A — AY\ipt)- 



IV. APPROACHING EQUILIBRIUM 

We now turn to our main topic, the approach to thermal equilibrium. Suppose we have 
an initial state = J2aCa\o:) such that the initial expectation value Aq = {iPq\A\iPq) of 
an observable A is far from its equilibrium value A = {A)t- What, then, is the behavior of 
At = (iptlAlipt) at later times? 

It is obvious that the answer depends on the details of the initial state. Without knowing 
them, we can only make a probabilistic analysis. There are two basic methods for doing 
so. One is to introduce a probability distribution for the initial state itself, and to average 
relevant quantities over it. The issue then becomes the justification of the procedure (e.g., 
maximum entropy) for choosing this particular distribution. The second method, which we 
will adopt, is to fix the initial state, and then study the values of interesting observables as 
functions of time. We treat the observation time as a uniformly distributed random variable, 
thus inducing a probability distribution for each observable. We then attempt to determine 
to what extent these probability distributions depend on the initial state. Ideally, there 
would be no dependence at all, thus rendering the choice of the initial state irrelevant. 

The moments of the probability distribution for At (which is induced by assuming a 
uniform probability distribution for t) are given by the infinite time averages 

lAtT =\im^ f^dt{AtY. (4.1) 

Again we note that we are using the infinite time average simply as a mathematical tool; in 
this context, the fact that astronomically long averaging times are necessary is not relevant. 
By using eq. (^.1|), and the nondegeneracy assumption discussed after eq. (|2.3|), we can 
express these moments as products of traces of powers of the matrix 

Aap = clA^pCp . (4.2) 

It will simplify the notation considerably if we first shift the operator A by a constant, so 
that the infinite time average of At is zero. This entails no loss of generality, and so from 
here on we will take 

A = Tri = 0. (4.3) 
The first few moments (A*)" can then be expressed as 
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{Atf = 2TiA\ 

'{A^ = 'OTi A^f + QTi A\ 

= 20 Tr Tr + 24 Tr I^ 

= 15(TrI2)3 + 90TrI2TrI^ + 40(Tri=^)2 + 120TrI^ (4.4) 

Determining these coefficients is a straightforward but tedious combinatoric problem; it 
involves counting the number of different ways that the time-dependent phases can cancel 
against each other, and so survive the infinite time average. Generalizing from eq. ( [4. 41) , 



the coefficient of a term of the form W^iTi A^^Y"^ in the expansion of {AtY is given by 
ral/dljpJw^D, where n = Y^ii^^iPi- 

We now wish to estimate the magnitude of Tr A™ for different values of m. Necessary 
inputs include the quantum uncertainty A in the total energy E, given by eq. ( p.5|) ; the 
energy bandwidth W of the matrix near energy given by 

\f(EM' ' ^ ' 

where f{E,uj) is defined in eq. (p.4|) ; the typical size a of the quantum/thermal fluctuations 
in A, given by 



a' ^ {A'Y = E (^')- = E \A^p\' ; (4.6) 

o a/3 



and the inverse participation ratio 



Z-ia \'-^a\ 



which counts the effective number of different energy levels that are present in the quantum 
state of the system. Also, X"^ can be regarded as the average value of Ic^p. Given X, it is 
helpful to define an effective level spacing between the participating energy eigenstates, 

SeS = J. (4.8) 

In general, 6es must be greater than or equal to the actual level spacing 6 ~ e^^E. 

Before proceeding to evaluate Tr A^, we must order the various energy scales. We expect 
that "typical" states of physical interest will have Ses ~ 5 and A ~ N~^/'^E ~ N^^'^T, since 
these are properties of a thermal density matrix. We also expect that W will be independent 
of N. To see why, we turn to the formula for At, eq. ( |3.1j ). From it, we can see that the time 
scale for the initial evolution of At is either h/W or h/A, whichever is larger. (Before this 
time, no relative phases have changed significantly.) For an observable A of interest, this 
time scale should be finite and nonzero in the thermodynamic limit, and hence independent 
of A^. Since we expect A ~ N^^'^T, it must be W that is independent of A^. We therefore 
conclude that 

4fT < < A < E (4.9) 
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is the regime of interest. 

With these considerations in place, we note that we can regard A as a banded random 
matrix of overall size {A/6es) x {^/^es) and bandwidth W/6es- Within the band of nonzero 
matrix elements, the magnitude of a typical entry is 



^typ ~ 



I-\W/5) 



^1/2 



a . 



(4.10) 



This comes from eq. ( [4.6| ), whose last equality demonstrates that a typical value of l^a/?] 
is {W/S)~^a'^; the extra factor of in eq. ( [4.10| ) comes from the Cq's in eq. ( ^^ . We can 
now estimate TrA™ for even m = 2k a.s 




(4.11) 



For odd m = 2k + 1, the expected value of TrA^'^"'"^ can be positive or negative; a root- 
mean-square estimate of its magnitude is 



TrA 



2fc+l 



A 



\ 2k+l 



typ 




2k+l / ^ \ {2fe+l)/2 




,2fc+l 



(4.12) 



However, this estimate should be regarded as less trustworthy than eq. ( [4.11| ). 

Eqs. ( ^.11|) and ( [4. 121 ) imply that the probability distribution for At has the form 



P{At) oc exp 




At 



At 



(4.13) 



where F^{x) and F_{x) are even and odd functions (respectively) whose Taylor expansions 
involve purely numerical coefficients. Eq. ( [4.13| ) can be verified by via a Feynman-diagram 
expansion for the moments; the quadratic term from provides the propagator, and the 
remaining terms (in both and F_) give the vertex coefficients. 

For At < {5/5e{iy^'^a, we can neglect F„ and all but the quadratic term in F+; we then 
have 



P{At) oc exp 



2 



A^' 



(4.14) 



where F^{x) = ^C,2x'^ + 0{x'^), and ^2 is a number of order one. Thus, for sufficiently 
small values of At, its probability distribution (as induced by a uniform distribution for t) 
is gaussian, and independent of the details of the initial state; only the energy width A of 
this state is relevant. Furthermore, we see again that the fluctuations of At about its inflnite 
time average (which is zero, by construction) are suppressed by a factor of (5/ A) 



1/2 



~ e" 



-S/2 



For At > (5/5efr) w, the nonuniversal behavior of the functions F± becomes relevant. 
Initial states that are "typical" according to most any reasonable criterion would all have 
5eff ~ S. In this case the nonuniversal corrections are important for At > a. 
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So far, our analysis has not addressed our original question: given the initial value Aq, 
what is the behavior of At at later times? To answer this question, we must compute the 
conditional probability P{At\AQ) to find the value At at time t, given the value Aq at time 
zero. By the usual rules of probability theory, this conditional probability can be expressed 
as 



P{At\Ac 



P{At,AQ 



(4.15) 



where P{At, Aq) is the joint probability for the observable A to have the expectation values 
At at time t and Aq at time zero. This joint probability is induced by assuming a uniform 
distribution for the initial time; hence the moments of P{At, Aq) are given by {At+t')'^ {At')"^ , 
where the time averaging is done with respect to t', with t held fixed. To compute these 
moments, we need an expansion analogous to eq. ( ^.4| ). Let us focus on the universal regime. 
In this case, the dominant terms in the expansion, as in eq. ([4.4|) , are those with the largest 
number of traces. Keeping only these terms results in a gaussian distribution. As always 
for a gaussian distribution, it is completely determined by its second moments. 



|2 



(A')' = E 1^-/31 



To streamline the notation a little we define a correlation function 

cit) = A^/'(A;y. 



(4.16) 



(4.17) 



Note that, by construction, C(0) = 1. 
These considerations lead to 



P{At, Aq) oc exp 



which has the correct second moments, 
conditional probability 



-I 









(At)' + (Aq)^ - 2C{t)AtAQ 



(4.18) 



[1 - C(t)]2a2 
Then, eqs. (|]T§, QO^I ), and (^J^) give us the 



P{At\Ao) (X exp 



2 



[At - C(t)Ao] 
[1 -C(t)]2a2 



(4.19) 



in the universal regime. Recalling that A/5 ~ e'^, eq. ( ^.19| ) shows us that it is overwhelm- 
ingly likely that we will find At = C{t)AQ. In the next section, we will see that the correlation 
function C(t) does not depend on the quantum state of the system. Thus, for all practical 
purposes. At is fully determined just by the initial value Aq; no other information about the 
state of the system is needed. 

For Aq > (6 /6es)^^'^a, the nonuniversal corrections can become important. We then 



expect a formula for P{At\AQ) that is similar to eq. ( |4.13| ). Because the large prefactor of 



A/5 ^ e should still be present, we still expect to get an effectively deterministic equation 
for At as a function of t and Aq, although it will no longer take the simple form At = C{t)AQ. 

Returning to the universal regime, we need to study the properties of C{t). This we do 
in the next section. 
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V. LINEAR RESPONSE 



Eqs. ( |4.1(j| ) and ( [4.17[ ) imply that the correlation function C{t) is real, even, and has a 
maximum magnitude of one. Also, it is quasiperiodic on the scale of the Heisenberg time 
th = 27Th/6 ~ e'^. However, this time scale is much too long to be of physical interest, and 
so we can justifiably ignore the quasiperiodicity of C(t). Then, also using eqs. (p.4|) and 
(ID, we get 

/ + 00 
du\f{E,u)\^e''^' . (5.1) 
-oo 

The last line shows that C{t) does not depend on the initial state, and that the bandwidth 
W of f{E,uj) sets the time scale for C(t). 

Eq. (|5.1|) can be compared with the results obtained through more standard methods. 
For example |2^, we can give A a nonzero expectation value at t = by supposing that. 



for t < 0, the system's hamiltonian was perturbed to H + XA, where A is a constant. For 
t < 0, we assume that the quantum state of the system is described by a thermal density 
matrix po ~ Q-iH+xA)/T ^ p^^, t > 0, this state is evolved forward in time with the original 
hamiltonian H, so that pt = e~*^*poe'^*; the time-dependent expectation value of A is then 
Af = Tr ptA. To leading order in A, this procedure results in A^ = CKubo(^)^O) where 

CKubo(t) oc ['^^ dp {AH{-ip)AH{t))T . (5.2) 
Jo 

Here Anit) = e*^*y4e~'^* is the operator A in the Heisenberg picture, and the angle brackets 
denote canonical thermal averaging as in eq. ( p.8| ). Eq. ( |5.2| ) can be written in terms of the 
matrix elements Aaf3 as 

^p^^ smh(c./2r) ^ (5.3) 

J-oo UJ 



This is not the same as eq. ( |5.1| ). Suppose, however, that the bandwidth W of f{E,uj) is 
smaller than the temperature T, and that the falloff of f{E,uj) for u ^ W is fast enough 
to make the integral in eq. (|5.3|) converge. Then the factor of sinh(co'/2T)/u; will be approx- 



imately constant over the relevant range of u, and so we will get 

CKuUt)=Cit) + OiWyT'). (5.4) 



In this case, eq. ( [4.19|) may be viewed as a verification of Onsager's postulate that a 



random fluctuation (in the value of some quantity) will dissipate in the same way as an 
initial value that is produced by an applied force. 

Let us now consider the circumstances under which the time evolution will be markovian. 
In the linear regime, and for t > 0, a markovian equation has the form {d/dt)At = —TAt, 
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where the parameter T must be real and positive. Given At = C(t)Ao, this imphes C(t) = 
exp(— r|t|). From eq. (|5?T|) , we see that this corresponds to \f{E,u!)\'^ oc + F^), and 

hence VF ~ F. However, l/iuj"^ + F^) does not fall off fast enough at large u for convergence 
of the integral in eq. (^.3|) , and so we expect some additional suppression of \ f{E, uj)\'^ when 
u is greater than the temperature T. If it happens that F is much less than T, then the 
time evolution will still be approximately markovian, but with "memory" effects on time 
scales less than h/T. This is consistent with other analyses P^ , |2^ , which find that the 
time evolution of expectation values is always non-markovian on time scales less than h/T, 
essentially because of the energy-time uncertainty principle. The overall conclusion is that 
markovian time evolution is associated with an isolated pole in \ f{E,u)\'^ at = ±iF, with 
F <^ T. However, this structure for \f{E,u!)\'^ is not required by any general principles, and 
so non-markovian behavior on all time scales is an open possibility. 



VI. DISCUSSION 



Our major result is eq. ( [4. 191 ), which gives the probability to find that an observable 
A has a quantum expectation value of At at time t, given that its expectation value is Aq 
at time zero (and assuming that A has been shifted, if necessary, so that the infinite time 
average of At is zero). This probability is computed for a particular quantum state l^pt), but 
with the zero of time treated as a uniformly distributed random variable. 

To understand the essential features of eq. ( [4.19| ), it is helpful to rewrite it more schemat- 
ically as 

P{At\Ao) oc exp {- O(e^) [At - C{t)Aof /a^} . (6.1) 

Here At = {^pt\A\^pt) is the time-dependent quantum expectation value of A; a} is the mean 
squared amplitude of the quantum fluctuations (which are also to be identified as thermal 
fluctuations) in A\ S is the thermodynamic entropy at energy E = {ipt\H\ipt)] and C{t) 
is the correlation function (normalized to one at t = 0) given in eq. ( ^.1|) , which, under 



favorable circumstances, is the same as the Kubo correlation function CKubo(''^) given in 
eq. ( ^.3|) . The range of validity of eq. ( |6.1|) is Aq < {5 / SesY^'^a., where 5 ~ E is the mean 
energy-level spacing near E, and 5efr > 5 is the mean level spacing for those eigenstates that 
participate significantly in \ipt)- Maximum participation (such as would be predicted by a 
thermal density matrix) corresponds to S^s ~ 5. 

Eq. ( |6.1| ) implies an effectively deterministic evolution equation for At, given only Aq as 
input: At = C{t)Ao. The probability that this equation is not satisfied is 0{e~^). Since 
C{t) is an even function of time, the effective evolution equation is invariant under time 
reversal. However, C{t) < C(0) for all t, and C{t) decays to zero if the bandwidth W, 
defined in eq. ([4.5|) , is finite. (There will be quasiperiodic resurgences of C{t) on the scale 
of the Heisenberg time th = 27Th/S ~ , but this is much too long to be relevant.) Thus 
a nonzero initial value Aq is overwhelmingly likely to evolve, both forwards and backwards 
in time, towards the equilibrium value of zero. This is of course entirely consistent with the 
heuristic picture of statistical mechanics originally proposed by Boltzmann. 

If the envelope function \f{E,u!)\'^, defined in eq. ( p.4|) , has an isolated pole at u = ±iT, 
then the time evolution will be approximately markovian, with C{t) ~ exp(— F|t|). This is 
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consistent with C{t) ~ CKubo(^) only if F <^ T, where T is the temperature corresponding 
to a thermodynamic energy of E. There will still be non-markovian "memory" effects on 
time scales less than h/T, which is consistent with the results of other analyses p8| , |2^ . 

To get eq. (|6.1j ) [or, more precisely, (|4.19| )], we had to make a fairly strong assumption 
about the quantum matrix elements of A, eq. ( p.4| ). However, this equation is well grounded 
in quantum chaos theory, and is likely to be at least approximately valid under a fairly wide 
range of circumstances. 

Naturally it would be useful to extend our results to the nonlinear regime, and naturally 
this is very much harder to do. Still we hope to return to this question in future work. A 
combination of our methods (which easily demonstrate the ubiquity of thermal equilibrium) 
with projection-operator methods (which easily generate equations for expectation values, 
assuming appropriately forced thermal density matrices as the starting point [^,^) might 
be a fruitful approach. 
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